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ABSTRACT 

Many stars exhibit stellar pulsations, favoring them for asteroseismic analyses. Interpreting the 
oscillations requires some knowledge of the oscillation mode geometry (spherical degree, radial and 
azimuthal orders). The 5 Scuti stars (1.5 - 2.5 M0) often show just one or few pulsation frequencies. 
Although this may promise a successful seismological analysis, we may not know enough about either 
the mode or the star to use the oscillation frequency to improve the determination of the stellar 
model, or probe the star's structure. For the observed frequencies to be used successfully as seismic 
probes of these objects, we need to concentrate on stars for which we can reduce the number of free 
parameters in the problem, such as binary systems or open clusters. We investigate how much our 
understanding of a 5 Scuti star is improved when it is in a detached eclipsing binary system instead 
of being a single field star. We use singular value decomposition to explore the precision we expect 
in stellar parameters (mass, age and chemical composition) for both cases. We examine how the 
parameter uncertainties propagate to the luminosity - effective temperature diagram and determine 
when the effort of obtaining a new measurement is justified. We show that for the single star, a correct 
identification of the oscillation mode is necessary to produce strong constraints on the stellar model 
properties, while for the binary system the observations without the pulsation mode provide the same 
or better constraints on the stellar parameters. In the latter case, the strong constraints provided by 
the binary system not only allow us to detect an incorrectly-identified oscillation mode, but we can 
also constrain the oscillation mode geometry by comparing the distribution of possible solutions with 
and without including the oscillation frequency as a constraint. 

Subject headings: asteroseismology — binaries: eclipsing — stars: fundamental parameters — stars: 
variables: delta Scuti — methods: analytical — methods: numerical 
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1. INTRODUCTION 

6 Scuti stars arc a class of pulsating objects located 
on the Hertzsprung-Russell (H-R) diagram on and near 
the main sequence (MS) where it intersects the classical 
instability strip (see Brcger 2000 for a review). They are 
1.5 - 2.5 Mq stars often pulsating in one dominant os- 
cillation mode or in many lower-amplitude modes. Orig- 
inally they were thought to be very interesting targets 
for an asteroseismic analysis because 1) their pulsations 
are easily detected with ground-based telescopes, and 2) 
we understand the structure of MS stars relatively well 
(e.g. Metcalfe et al. 2009). Consequently, asteroseismol- 
ogy of these stars was thought to have the potential to 
probe the details of the interior, such as energy trans- 
port mechanisms and convective core overshoot, as well 
as less well-understood phenomena such as rapid r otation 
(|Featherstone et al.l[2007l : iMacGregor et al.ll2007( ). 

Unfortunately, successful asteroseismic analyses are in- 
frequent for several reasons. 1) Observationally, the 
labeling of each measured frequency with its geomet- 
ric characteristics (mode identification) is not trivial. 
For example, rapid rotation causes each of the mode 
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degrees £ to split into {2i + 1) components w ith dif- 
ferent azimuthal order (m) (jGoupil et al.l 120051 ). Also 
some o f the stars are suffici ently evolved to show mixed 
modes (jMetcalfe et al.|[2010() , thus complicating the anal- 
ysis. 2) Theoretically, it is difficult to find a unique 
model to match the set of observed frequencies, even 
with additional observational constraints. With the 
fundamental stellar properties poorly constrained, use 
of the mode as a "seismic probe" is severely limited. 
For example, the multi-periodic 5 Scuti star FG Vir, 
has been the target of many observational campaigns, 
allowing the detec tion of more than 75 frequencies 
Breger et al.l 120051) . and the identif i cation of 12 mode s 



Daszvhska-Daszkiewicz et al.l 120051 : IZima etlill |2006|) . 



Even with these observational constraints, the stellar pa- 
rameters s uch as mass and helium mass fraction remain 
uncertain (iGuzik fc BradlevI 119951: IViskum et al.l 119981: 



iTempleton et al.l 120011: iKirbivik et all I2005D . A simi- 
lar case is XX Pyx fsee iHandleii 120021 for a review). 
iHandled (I2002D noted that even searching among more 
than 40,000 models in a three -dimensional parameter 
space, iPamvatnvkh et afl (jl998[ ) could not find a model 
that matched the observations. 

The lack of successful asteroseismic analyses has 
encouraged authors to reconsider their approach. 
On the one hand, much progress has been made 
in observational mode identification, while on the 
other, authors are beginning to study these stars 
in syste ms where the number of free para meters is 
reduc e d (iLampens fc Boflin 2000l: lAerts &: Harmanc3 
120041: iMaceronil 120061: ICosta et al.l I2007D . For ex- 
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ample, open cl usters provide constraint s on the age 
and metallicity ()Fox Machado et all l2006l : ISaesen et alJ 
|2010() . while multiple systems allow precise determina- 
tions of the mass and radius (lEscola-Sirisi et alJ 120051 : 
iHekker et alll2010t iMaceroni et alJl2010( ). The ideal pul- 
sating star wo uld be part of an eclipsing binary within an 
open cluster (lArentoft et alJ 120071: ICreevev et all l200l 



IStello et al.ll2010l:lTalamantes et al.ll2010[ ). 

Using asteroseismology to probe the interior of a star 
requires the global properties to be known quite well. 
For example, the mass should be known to 1-2%. Fortu- 
nately, the measurable quantities from a binary system 
can provide strong constraints on the properties of the 
component stars. If the binary is an eclipsing and double- 
lined spectroscopic system (SB2), the absolute values of 
the masses and radii can be d etermined with 1-2% p re- 
cision e.g. iRibas et all (|1999[) : iLastennet et all (|2000[ ). 

Many new detached eclipsing binaries are being 
discovered by several satellites dedi cated to photo- 
mctric monitoring of sta rs fMO ST iMatthewsl I2004L 



CoRoT iBaglin et all [20021 Kevler iBorucki et al 
well as numerous ground-bas e d surveys 



as 



(iHenryl 

19981 iBrown fc Charbonneaul 120001 : iBakos et al.l 120021 : 
PoUacco et al.ll2006D . Some of these sys tems have com 
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poncnts that exhibit stellar oscillations ([Maceroni et al.l 
[2010; Hckkcr ct al. 201(1). We wiU soon be faced with 
the problem of choosing which systems to study. 

In this paper, we use singular value decomposition 
(SVD) to study the information content of stellar sys- 
tems. This theoretical investigation quantifies the in- 
crease in astrophysical information realized by studying 
pulsating stars in detached eclipsing SB2 binaries (bi- 
nary) instead of single pulsating field stars (single star), 
especially for those cases where mode identification is 
difficult. We assume that in a detached system, the stel- 
lar structure is unaltered by the binarity. The particular 
questions we address are: 

- What is the role of each observable for a single star, 
and how do these roles change if the star is observed 
in a binary system? 

- With the additional constraints for the binary sys- 
tem, how does the error box in the H-R diagram 
compare to that of a single star? If the mode geom- 
etry is successfully identified through observational 
methods, how does this error box change? 

- Given the typical observational errors, what are the 
expected uncertainties in the model parameters for 
the single star and the binary system? 

- Are the constraints on the stellar parameters suf- 
ficient to distinguish between possible solutions 
when a mode is incorrectly identified? 

In Section [2] we introduce the definitions, models, ob- 
servations and method necessary to understand the rest 
of this work. In Section [3] we study the roles that the 
observations play for determining the stellar model solu- 
tion for both a single star and a binary stellar system, 
and we highlight which observations are most informa- 
tive. In Section |3| we then discuss the parameter un- 
certainties, the error ellipses, and the error box in the 



H-R diagram, which show more specifically which obser- 
vations are capable or incapable of reducing these uncer- 
tainties. Finally in Section [S] we use simulations to show 
that an incorrectly-identified oscillation mode can be de- 
tected in the binary system, but not in the single star. 
We also show that the binary observations alone can con- 
strain and even identify the mode geometry by studying 
the distribution of model solutions without observational 
constraints on the mode. 

2. MODELS, OBSERVATIONS AND METHODS 

Our method follows that of iBrown et al.l (|1994[ ) who 
used SVD as a diagnostic tool to investigate how use- 
ful oscillation frequencies arc as constraints on stel- 
lar parameters in a visual binary system, using a Cen 
A as an exa r nple. Th is work was f ollowe d up in 
ICreevev et all (|2006D and ICreevev et all (I2007D an d the 



techniques are also present ed in [Press et alT l 19921 ) and 
iMiglio fc MontalbanI (|2005[ ). We refer the reader to these 
papers for details on the method. 

2.1. The distinction between parameters and observahles 

We define a parameter as an input characteristic to 
a stellar model, for example, mass and age. These are 
the quantities that we wish to determine. An observable 
is an output quantity from a stellar model given a set 
of parameters e.g. a radius, metallicity or a photometric 
magnitude V. We compare the theoretical observables to 
the real observations in order to retrieve the parameters. 
For example, we can measure the effective temperature 
of a star (observation) and its error, and use these to 
retrieve the mass and age (parameters) for a stellar evo- 
lution model by comparing the observed effective tem- 
perature to the theoretical/model effective temperature 
(observable) . 

2.2. Stellar parameters and models 

We describe a single 5 Scuti star by a set of parame- 
ters P. These are mass M, age r, rotation (we use ro- 
tational velocity w), initial hydrogen X (or helium Y) 
and heavy metal Z mass fraction where X -\-Y + Z = \, 
and mixing-length parameter a, where applicable. Fig- 
ure [T] shows the derivatives of the two most sensitive 
observables (the effective temperature and an oscilla- 
tion mode) with respective to a for a range of stel- 
lar masses. For values above 1.7 M0, we can clearly 
see that the considered observables are insensitive to 
a, and can be ig nored. The Aarhus STellar E volution 
Code (ASTEC) ([Christensen-Dalsgaardl [20080 is used 
to calculate the stellar evolution models . The ASTEC 
code uses the equation of state of Egglet on et al.l (|1973f ). 
without Coulomb corre ctions, and the OPAL opacities 
l|Igelsias fc Rogerj|1996f ). supplemented by Kurucz opac- 
ities at lo w temperatures. The nuclear reaction rates 
come from lBahcall fc Pinsonneaultl (| 19921), convection is 
descri bed by the mixing- length theory of iBohm- Vitensi 
(|1958[ ). convective core overshooting is included with 
aov set to 0.3, and diffusion effects are ignored. This 
code uses the stellar parameters described above as the 
input ingredients and returns a set of global stellar 
properties such as radius R and effective temperature 
Tcff (observables), and the interior profiles of the stel- 
lar mass, density and pressu re. Oscillation frequencies 
are calculated using MagRot ([Cough fc Thompson|[l990l : 
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Fig. 1. — The dependence of the most sensitive observables 
(effective temperature and an oscillation mode) on the convective 
envelope mixing-length parameter for a range of stellar masses. 



TABLE 1 
System Parameters 

Parameter Value [Pj) Parameter Value (P,) 



Ma 


1.80 Mq 




100.0 km s-i 


Mb 


1.70 Mq 




60.0 km s^i 


X 


0.700 


Z 


0.035 


T 


0.700 Gyr 


d 


200 pc 


a 


0.150 AU 


i 


85.6° 


e 


0.0 


7 


0.0 km s~i 


(a; 


0.0) 







iBurke fc ThompsonI |2006[ ). The distance d is also in- 
cluded as a stellar parameter, and coupling this with R, 
M, Tefi and the metallicity [M/H], allow s us to calcu- 
late m agnitudes using model atmospheres (jLeieune et alj 
[1991 . 

For a binary system, the additional model parameters 
are the system properties: orbital semi-major axis a, or- 
bital eccentricity e, longitude of periastron w, systemic 
velocity 7, and orbital inclination i (note that i will al- 
ways denote inclination unless otherwise specified). Both 
stellar components of a binary system share the param- 
eters T, X and Z, so the individual stars differ only in 
M and v. We shall use the subscripts 'A' and 'B' to de- 
note the components of the binary. The parameters of 
the model system are given in Table [U and while these 
are not based on any particular binary system, there are 
known systems whose masses approximate these, such as 
HD 172189, HD 26591, HD 42083. We also note that 
the system stars pursue non-synchronised rotation, com- 
patible with our earlier assumption that the stars evolve 
" individually" . 

2.3. System observables and observations 

The observables are the measurable quantities of the 
system. These include things such as R (from interfer- 
ometry for example), T^s, [M/H], gravity logg, and par- 
allax TT for a single star. For a binary system, additional 
observables include effective temperature ratio Tba (= 
Tb/Ta), relative radii Rba, semi-amplitudes of radial 
velocity curves Ka and Kb, and orbital period IT. We 
note that the radial velocity or light curve measurements 
are also considered observables, however, in this work the 



more intuitive derived values (Tba Ma sin"^ z, Mb sin'^ i, 

oj etc.) are used. 

Using model atmospheres (jLeieune et al.|[l997l ) with d 
and i?, we calculate the photometric magnitudes (r, t, 
and z) using the SDSS filter system l|York et al.l l200C/). 
These atmospheres take reddening into account. We sub- 
sequently calculate the colors from these magnitudes, al- 
though we note that distance d is not necessary to ob- 
tain approximate color indices. The colors we include in 
our analysis arc (r — i) and (i — z) and the magnitude 
is r. For the binary, the photometric observables are 
calculated from blended spectra. Two or more blended 
colors yield estimates of the individual and relative quan- 
tities of the component properties. The c olors of the 
ind ividual comp onents can b e disentangled (jKochl [19601 : 
iSemcn iuk 2001i: [Creevev et'al . 2005). However, obtain- 
ing component colors with the same precision as the sys- 
tem can only be achieved with high quality multi-color 
light curves. 

From the light curve of an eclipsing binary system we 
obtain Rba, Tba, i and H. We use these derived values 
from the modeling of the binary system for simplicity 
and to understand the errors in terms of what is quoted 
in the literature. The derived observables are valid, as 
long as we correctly propagate the uncertainties from the 
light curve data. 

Spectroscopy provides the observables va sin i and 
VBSmi (projected rotational velocity) and the atmo- 
spheric parameters logg. Toff and [M/H]. An observed 
spectrum for a binary system is the sum of the spectra 
from the two stellar components. Many observational 
techniques exist to disentangle spectra e.g. iHadraval 
(|l995l,[2009l). and these techniques have been successful, 



Creevev et al.l ()2009l ) for a specific case using various 



see 

techniques to determine the individual component Toff 
for HD 172189. We assume that we can successfully dis- 
entangle the spectra and obtain effective temperatures 
for both stars while using the photometric Tba as a con- 
straint. 

A time series of spectroscopic measurements yields a 
time series of radial velocities, and modeling these data 
yields Masiv? i, M^sin'^z, asini, e, lo and 7. Coupling 
these values with i and Rba (from photometry) provides 
the absolute values of the radius. 

A time series of radial velocity or photometric mea- 
surements provides the frequcncy(ies) u of the modes 
that are excited in the star. However, the frequencies 
are only useful if wc know the oscillation mode geometry 
(the degree I, the azimuthal order m and the radial or- 
der n of the wave) . This so-called " mode- identification" 
can be achieved by either studying the time series of 
spectroscopic line-profiles of the star (iKennellv et al 



199C; 



20031: 



le-pronles or tne star llKen 
Horner et al.lll996t iBaloMi 120001: [Bri^et fc Aerts 
Zimall2006f) . bv using multi-color photometry and 



comparing the oscillation amplitudes an d the phases of 
the frequencies in various filter bands (iGarridd l2000l : 
iDupret et"al[l2003t iDaszvns ka-Daszkiewicz et al.l |2007D, 
or by using t he screening effects of the eclipses in eclips- 
ing binaries (jBiro fc Nusp]||2005t [Camarova et al.| [2005l. 
Many 5 Scuti pulsators have only one dominant oscilla- 
tion frequency and so the assumption is that we measure 
at least one precise frequency. 

The values of the errors (c.f. Table[2]) arc chosen by con- 
sidering typical observations for such systems, but also 



4 



Creevey et al. 



TABLE 2 

Observables and Errors 



Observable {Bi) Error (e^) 

Figs.[2T3l5l Figs. I4l8l Figs. [Tim 



Ra (Re) 


1.95 


0.02 


0.01 


0.04 


Rb (R©) 


1.81 


0.04 


0.01 


0.04 


Tba 


0.97 


0.05 


0.05 


0.05 


i (°) 


85.6 


1.0 


1.0 


0.3 




6965 


100 


100 


100 


va sini (km s~^) 


99.7 


2.5 


2.5 


2.5 


Vb sini (km s~^) 


59.8 


2.5 


2.5 


2.5 


n (yrs) 


3.1053C-6 


lc-6 


le-6 


le-6 


Ma sin^ i (Mq ) 


1.78 


0.06 


0.10 


0.04 


Mfl sin^ i (Mq) 


1.69 


0.05 


0.10 


0.03 


[M/H] (dex) 


0.31 


0.05 


0.05 


0.05 


logs (dex) 


4.1 


0.5 


0.5 


0.5 


TT (mas) 


5.0 


0.5 


0.5 


0.5 


r (mag) 


8.16 


0.02 


0.021' 


0.01 


(r - i) (mag) 


0.178 


0.003 


0.005t 


0.001 


(j - z) (mag) 


0.009 


0.003 


0.005t 


0.001 


(/iHz) 


1108.9 


1.3 


1.3tt 


1.3 



t 0.5 mag for Fig. [4] top panel. It ie4 ^Hz for the black ellipse in 
Fig. [8] 



to study some limiting cases, for example when the er- 
rors reach the limits of current techniques, or if the stars 
are bright (closer in distance). We assume that the light 
curves cover a few orbital periods. Spectroscopic data, 
however, are more time-consuming and the demand for 
observation time is high, especially for larger telescopes. 
Thus, one may hope that spectroscopy during one full 
orbital period can be obtained, with some data points 
during subsequent orbits. Taking into account that one 
of the stars is pulsating, this will impact the analysis of 
both the photometric and spectroscopic data. In partic- 
ular the depths of the eclipses will need to be observed 
several times to remove the effects of the modulations due 
to pulsation, and this influences the derivation of Tba, 
Hba, and 11, for example. Spectroscopy will show clear 
line-profile variations, and deviations in the expected ra- 
dial velocities due to pulsation. Both of these will influ- 
ence the determination of the orbital radial velocities to 
some extent, and these will in turn affect the errors in 
M, R and 11. Tabled first column, lists the values of the 
observables calculated from the stellar model defined by 
the parameters in Tabled] The subsequent columns show 
the corresponding errors used for this analysis. We have 
considered errors (e) in mass-related and radius observ- 
ables of the order of 1-3%. The 3% error should account 
for the effects of pulsation on the data. Errors of 100 K, 
0.05 dex and 0.3 - 1.0° for Teg, Tba, and i are values 
typically found in the recent literature ("Ribas et al.ll2002l: 
Southworth et al. 2004; Hilditch et al. 2005; Sousa et al.l 
20061: iBruntt et"arH2010D . We use 2.5 km s "^ as the error 
on th e projected rotational velocity (e.g. ICreevev et al.l 
I2009| ). Reducing the precision on these values influences 
the determination of the rotational velocity of the star. 
Unfortunately including rapid rotation is a complex mat- 
ter for stellar models, and most stellar codes (including 
those used in this work) consider only solid-body rotation 
which has minimal effect on the non-seismic observables. 



TABLE 3 

Observable (OS) and parameter (P) sets used in this 

ANALYSIS. AS DENOTES ADDITIONAL OBSERVABLES AND THE 
subscripts 'S' AND 'EB' REFER TO SINGLE STAR AND ECLIPSING 
BINARY SYSTEM. 



OSs = {fl,T3ff,log(;,i;sini,[M/H],7r} 

OSeb = {RA,RB,TBA,i,TefiA,VAsini,VBsmi, 

n, Ma sin^ i, Mb sin^ i, [M/H], vr} 
ASl = {u} 
AS2 = {r,r — i,i — z} 
ASS = {r,r — i,i — z,u} 

Ps = {M,T,X,Z,v,i,d} 

Feb = {MA,MB,T,X,Z,VA,VB,->;d, [a,e,uj,'y]} 



1.0 
0.8 
0.6 
0.4 



HP 0.8 

m I 1 

0.6 I I 

0.4 I I 

0.2 

0.0 I II II II II II II 

R Teff logg vsini[M/H] tt r r — i i — z 

Fig. 2. — Significance of observables sets without colors, magni- 
tudes and the oscillation mode (OSg — top) and including these 
data (OSs+AS3 — bottom) for a single star. 



The precision in the photometric values is typical for 
well-determined literature values. We chose some lim- 
iting values (e < 0.003) because at these precisions the 
colors begin to have an important impact for the deter- 
mination of the parameters. The frequency error of 1.3 
fiHz corresponds to the frequency resolution of data from 
a one- week observational campaign. We increase the val- 
ues of these errors for some calculations (note in Table [5]) 
to eliminate their influence. 

An analysis of the influence that each observation has 
(see next section) for different sets and errors, and an 
investigation of the parameter uncertainties led us to 
choose four competitive sets of observables to study, and 
these are given in the top part of Table [3l OS is the 
base observable set used, and ASl, AS2 and ASS are 
the additional sets. The subscripts 'S' and 'EE' denote 
'single star' and 'eclipsing binary' respectively. For the 
remainder of this work we discuss our results in terms 
of each of these sets and we remove the bold text: OS, 
OS-^-ASl, OS4-AS2 and 0S-^AS3. We include the sub- 
scripts 'A' and 'B' on the observables and parameters to 
denote each component in the binary system when neces- 
sary, and the magnitudes and colors also have a subscript 
'EB' to emphasize that these are binary system (blended) 
observables. 

We do not include a sin i as an observable because it 
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does note provide independent information. We can cal- 
culate it from three of the following: Ma sin^ i, Mb sin^ 
q = Ka/Kb 
Law. 



Mb /Ma, and 11 using Kepler's Third 



2.4. Singular value decomposition method 

SVD is the decomposition of any M x N matrix D into 
3 components: U, V and W given by D = UWV'^. 
V"^ is the transpose of V which is an x orthogonal 
matrix that contains the input basis vectors for D or 
the vectors associated with the parameter space. U is 
an Af X iV orthogonal matrix that contains the output 
basis vectors for D, or the vectors associated with the 
observable space. W is a diagonal matrix that contains 
the singular values of D. 

The key element to this work is the description of the 
matrix D, which we call the design matrix and each el- 
ement is a partial derivative of each of the observable 
with respect to each of the parameters of the system, 
taking into account the measurement errors for each of 
the obscrvablcs: 

Dk, = ^/e,. (1) 

Here are each of the k = 1,2, ...M observables of the 
system with expected errors e^, and Pj are each of the 
j = 1, 2, free parameters of the system (see section 
l2.2l for discussion on the observables and the parameters). 

By writing the design matrix with the measurement 
errors taken into account, we provide a quantitative de- 
scription of the information content of each of the ob- 
servables for determining the stellar parameters and their 
uncertainties. 

Starting from an initial close guess of the solu- 
tion Pq, SVD can be used as an inversion tech- 
nique to obtain the true solution Pr of the system. 
This is done by calculating a set of parameter correc- 
tions (SP_that minimizes some goodness-of-fit function: 
5P = VW^-'-U'^(5B, where (5B are the differences be- 
tween the set of actual observations O and the calculated 
observables Bq given the initial parameters Pq. W is a 
modification of the matrix W such that the inverses of all 
values below a certain threshold are set to 0. The formal 
errors comprise the sum of all of the Yq/wq, where each 
Wq/wq describes the direction and magnitude to change 
a parameter so that the true solution Pr and formal un- 
certainties can be given by 



VW-^U'^^B I ± 



Vi , V: 



± 



± ... ± 



V 



N 



Wi 



W2 



WN 



(2) 

The covarianee matrix C consequently assumes a neat 
and compact form: 



N 



= E 

9=1 



VjqVlq 



(3) 



and the square roots of the diagonal elements of the co- 
variance matrix are the theoretical parameter uncertain- 
ties. Note that e is reserved for observational error and 
<7 for parameter uncertainty: 



N 

E 

9=1 



J9 

Wn 



(4) 



Another useful property of SVD is the significance S of 
an observable. This can be quantified as a measure of the 
extent that a 1-e change in shifts the inferred param- 
eters towards the l-cr error ellipsoid in parameter space. 
In this way, S quantifies the impact that an observable 
has for the determination of the parameter solution: 

1/2 




(5) 



Because of the orthonormality of the decomposition ma- 
trices, the value of S varies between and 1. A low Sk 
implies that the observable B^ has relatively less influ- 
ence on the solution, and that a change in the measure- 
ment will have little or no impact. A high value of Sk 
implies that this observable is important for the solution 
and any change in the measurement will force a corre- 
sponding change in the solution. 

Finally, the matrices U and V provide information 
about the role each of the observables plays in deter- 
mining the stellar parameters. Each column vector of U 
[Uj) is related uniquely to each column vector of V (V^) 
and its importance in the solution is given by the corre- 
sponding singular value Wj (cf. Fig Sj) . In section 13.1.21 
we elaborate on this discussion. 

3. BINARY SYSTEM VERSUS SINGLE STAR 
CONSTRAINTS 

In this section we investigate the roles that each ob- 
servable plays for determining each parameter. We study 
how these roles change when we include/exclude certain 
observables and when we consider different values of the 
observational errors. 

3.1. Single star observables 
3.1.1. Significance 

Figure [2] shows the significance (Eq. [5]) of each of the 
observables for a single star for OS5 (top panel) and 
OS5+AS3 (lower panel), using the measurement errors 
from the second column in Tabled Because there is no 
information about i (inclination), the figure illustrates 
that vsini is not an effective constraint {S = 0). Al- 
though, the value of vsmi imposes a lower limit on v, 
we can not determine its uncertainty, log g is also a weak 
constraint, because both R and T^s provide similar in- 
formation, but of better quality. If either of these were 
not available, then \ogg would have a higher significance. 
5 = 1 implies that R, Teff , [M/H], and tt are all necessary 
observations to constrain P5. By including the photo- 
metric information and the mode (lower panel), the most 
noticeable change is the reduction in S'(rcff), whose infor- 
mation is superseded mainly by r and {r — i). S{v) = 1 
indicates that the mode has a strong influence on the 
determination of the solution. 

Figure [3] shows how S changes for each observable as 
some measurement errors arc varied. The left panel 
shows S{r - i), S{R) and S{T^s) varied for 

0S5-hAS2. Similar results are found for OS5-fAS3. 
S[r — i) decreases as its error increases to values of about 
0.004 mag, and, at this value, S{R) begins to supersede 
S{r — i). S{Tcs) also increases as the color measurement 
is more poorly determined, although it never increases to 
more than 0.45. The largest increase (or change) is seen 
only as e^-i reaches mmag precisions. 
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Fig. 3. — Change in observable significance S as the precision in 
photometric information changes for a single star. The photometric 
observables are denoted by crosses on their lines. Left; Similar 
results are found for sets OSs+AS2 and OSg+ASS, as the precision 
is decreased in e^-i. Right: Continuous and dashed lines denote 
OSs+AS2 and OSs+AS3, respectively, as the precision in is 
decreased. 



The right panel of Figure |3] shows the change in 
5'([M/H]) and S{i — z) as ei-z increases in value. In this 
case, no other observable is affected, {i — z) has a strong 
influence for determining the parameter solution only at 
precisions of < 0.002 mag; for (r — i) the critical value is 
^0.004 mag (left panel). The difference between the con- 
tinuous and the dashed lines is the inclusion of seismic 
data (ASS, dashed). Without seismic data, both [M/H] 
and {i — z) are important observables, while including 
seismic data, [M/H] has little impact when [i — z) is im- 
portant, and vice versa. It is primarily the parameter Z 
that these observables arc responsible for determining. 

3.1.2. Transformation matrices 

The decomposition matrices of SVD show which ob- 
servables determine each of the parameters of the system. 
These are sets of linear vectors where each vector in U 
is related uniquely to each vector in V, and the relative 
importance of each vector is given by the correspond- 
ing singular value Wj in W. In Figure 01 we show the 
transpose of these decomposition matrices for the single 



star system, with the singular values represented in the 
right-most figure of both panels. These figures take into 
account the errors given in the third column of Table [U 
except for the upper panel where the values for the colors 
have been increased to 0.5 mag and the frequency error 
to le4 //Hz to eliminate their influences. 

In FigurelHtop panel, R, Teg, and [M/H] appear in the 
top three observable vectors indicating their importance 
for determining the combination of parameters, shown in 
vectors Vi, V2 and V3 (mass, age, chemical composition). 
The next important observable (in fjg — here wq > W4) 
is TT, and this is uniquely responsible for determining the 
distance (in Vg). Likewise, only the observable wsini de- 
termines V, but poorly, and finally the observable with 
least information is \ogg {U3 with smallest 1^3), which 
contributes weakly to determining mass, age and chemi- 
cal composition. 

In the bottom panel of Figure |4l we have improved 
the errors in the photometric observables and the oscil- 
lation mode (third column in Table [2]). These new ob- 
servables now play a dominant role for determining the 
parameters, as can be seen from their positions in the 
top vectors, as well as an increase by a factor of 10 in 
the top singular values. Mass, age and chemical compo- 
sition are much better constrained, mostly appearing in 
the top two panels. The contributions from Tcn and log g 
have diminished considerably, and the distance is now no 
longer determined uniquely by tt. It also depends weakly 
on the photometric colors and magnitude (V5 and U5). 
As in the top panel, v sin i is responsible for determin- 
ing V (Uq and Vg). By omitting the oscillation frequency 
from the set of observables. R becomes the observable 
mostly responsible for determining M. 

3.2. Binary system observables: Significance 

The top panel of Figure [5] illustrates S for OS^s 
(binary system) using the errors from the second col- 
umn of Table [2) This figures shows that the most 
important quantities are Ra, Rb, Tcsa, H, [M/H], 
Ma sin"^ i, Mb sin^ i and tt. The lower panel shows S for 
OSsB-f ASS. Note how each of the S values changes: in- 
cluding the magnitude, colors and oscillation frequency 
reduces the importance of Ra, Tcsa, [M/H] and Tba, 
while ever so slightly reducing 5(7r). We note that i ap- 
pears deficient in information content. An accurate and 
precise determination of i is necessary to correctly scale 
the observed values of Ma sin'^ i and Mb sin^ i. In this 
case e(z) = 1.0°, however, increasing its precision to < 
0.3° causes S{i) to increase rapidly. 

Figure [6] left panel shows S{Tcsa) and S{r — i)EB as 
^r-iEB changes for OSbb-I-ASS. As er-iSB decreases in 
precision, S{r — i)EB decreases from 1.0 to 0.6, while an 
increase in S{Tcsa) from 0.6 to 0.8 is also seen, perhaps 
implying that (r — i)EB is a temperature indicator. This 
is similar to the single star case but in the binary sys- 
tem Toff remains relatively more important because the 
blended (r — i) eb provides less information (see Fig. [3] 
left diagram). 

The right half of Figure [H] shows S{Ra) as its pre- 
cision is decreased for OSeb and OSes-fASl. In the 
latter case, including an identified mode causes S{Ra) 
to rapidly decrease in value, as it becomes more poorly 
determined. For some observables (Tba, Ma. b sin^ i) S 
increases slightly during this change. At more than about 
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Fig. 4. — Decomposition matrices of a single star model using the errors given in Table [2] third column. The left, center and right panels 
are the matrices , V""^ and W respectively. The difference between the panels is the addition of well-constrained photometric and seismic 
data in the lower panels. These values are set to 0.5 mag and le4 /iHz respectively in the upper panel. 
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2% error. Ra has little role to play. On the other hand, 
for OSeb, Ra remams relatively important at all consid- 
ered errors, simply because there is no other observable 
in OSeb that provides the same type of information. 

4. DETERMINATION OF THE FUNDAMENTAL 
PARAMETERS OF THE PULSATING STAR 

The observables that showed most sensitivity to a 
change in their precisions for different sets for either a 
single star or a binary system are R, Tgff, the photo- 
metric colors and an identified mode. In this section we 
study the uncertainties (Eq. 2]) in the global parameters 
of the (primary) pulsating star as these observables are 



improved/included for both a single star and a binary 
system. 

4.1. Parameter uncertainties 

To successfully use a single identified mode as a "seis- 
mic probe" , we need to know the fundamental stellar 
parameters M , t, X, Z and v of the pulsating star with 
good precision. Here we discuss only the uncertainties in 
the first four, because the determination of v is indepen- 
dent of the determination of the others (the covariance 
off-diagonal elements ^ 0) and independent of most of 
the observable errors, except i and w^sini. All of the 
other parameters of the system {Mb sin^ i, e, etc.) are 
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Fig. 7. — Theoretical uncertainties in the stellar parameters as the precision in en (left) and ei_2 (right) increases, for a single star 
(dashed lines) and a binary system (solid lines). The lines with dots represent the sets of observables with seismic information. The left 
panel shows the results without photometric constraints. The y-axis for each parameter has the same scale for comparing panels. 



included for deriving the uncertainties, due to correla- 
tions that must be accounted for. However, we do not 
discuss these uncertainties because we are only interested 
in describing the pulsating star. As we are dealing with 
the primary (A) pulsating star only, we drop the sub- 
script 'A' on all of the quantities. 

Figure [7] shows the theoretical uncertainties in M 
(top), T (second), X (third) and Z (lower panels) as we 
increase the precision in R (left panels) and both of the 



photometric colors (right panels). The latter is denoted 
by ei^z and for the binary it implies the blended sys- 
tem's colors. The dashed and solid lines represent the 
single and binary system respectively, and the observ- 
able set with the identified mode has dots overplotted on 
the lines. The results for OS (both OS5 and OSeb) and 
OS-l-ASl are represented on the left panels, and those for 
0S-I-AS2 and 0S+AS3 (including photometric informa- 
tion) on the right. Both the left and right panels show 
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the same y-axes scales for comparison. 

It is very clear that the oscillation mode plays a very 
different role in the single star and the binary system. 
Excluding it from the single star constraints, without 
photometry, yields uncertainties > 100% and indepen- 
dent of en', the constraints are inadequate to determine 
the stellar model. Including the color constraints also 
yields large uncertainties (> 30%) although not entirely 
independent of ti-z- When the mode is included, the 
uncertainties improve rapidly as the errors in both types 
of observables improve. In fact, cr{M) becomes strictly 
dependent on eji as the error in the radius improves to 
< 2 - 3%. 

There is little difference seen in the parameter uncer- 
tainties for the binary system when the oscillation mode 
is included/excluded as a constraint. The binary sys- 
tem without any seismic data provides similar or better 
constraints on all of the parameters than the single star 
with an identified mode, except in some exceptional cases 
(see below) . The fact that including a mode with the bi- 
nary system constraints does not lead to improvements 
in the pulsating star's parameter uncertainties indicates 
that the binary constraints alone may provide sufficient 
constraints on the model to use the identified mode in a 
different way. 

For the binary system, we also tested the effect of im- 
proving the precision in Ma sin'^ i and Mb sin'^ i. Setting 
both of these values to 0.01, 0.05, and 0.10 M© (0.5%, 2% 
and 5%) has little effect: for = 2%, ct{M) ^ 0.3, 1.0, 
and 1.5% respectively with and without the identified 
mode, and for eu > 4%, a{M) levels out to 0.3, 1.5, and 
2.5%. 

We note that for both a{Z) and CT(r) for 0S+AS3, in 
some cases the single star provides slightly better con- 
straints than the binary. This happens only for very 
precisely measured (unblended) colors, while including 
an oscillation mode as a constraint. Finally, one should 
also note that a{X) > 10% is not a good constraint. The 
absolute value of the parameter is 0.700, and a 0.070 er- 
ror on this value produces no meaningful constraint on 
X. 

4.2. Parameter correlations 

The V matrix from the SV decomposition describes 
the correlations among the parameters, while the sin- 
gular values scale each of the vectors Vj to produce the 
n-dimcnsional error ellipse axes. In Figure[S]wc show pro- 
jected contours of the l-cr error box for r versus M, and Z 
versus r using the vectors that yield the largest ellipses. 
The black/grey ellipses correspond to the sets of errors 
in the third column of Table [2] for OS£;b+AS2/AS3. In- 
cluding the oscillation mode does not reduce the error 
ellipse, and this reinforces the possibility of using this 
extra information to learn something else about the star. 

4.3. H-R error box 

Using equation ^ we calculate the theoretical uncer- 
tainties in the model effective temperature and model lu- 
minosity. Figure ^ shows the approximate theoretical er- 
ror boxes for these quantities for a single star (OSg-l-ASl, 
dashed lines) and a binary system {OSeb solid lines). 
For the single star an identified mode is included because 
Fig. [7] and Sect. 14.1] showed that the parameters are not 



constrained if the identified mode is not included, while 
for the binary system no seismic data are included. 

The error box shrinks significantly in the luminosity 
axis while reducing the error in the radius for the single 
star (left panel). cr(i*) also changes slightly with ej-ott 
(right panel), but its value is determined primarily by the 
error in the radius (2%). The left panel shows that in- 
terpolating between 1% and 3% for cr produces a a■{L^,) 
= 0.5 L© for en = 2%, as shown in the right panel. 

As expected, the error box shrinks in the model effec- 
tive temperature axis due to the reduced errors in the 
observational T^g (right panel). The observational et^jj 
of 200, 100, and 50 K produce model uncertainties cr(reff) 
of 250, 110, and 50 K. We do not expect these values to 
be reproduced exactly because the observed Tcff has a 
measurement error associated with it, and for large val- 
ues of the error ('-^200K), other observables can domi- 
nate the determination of the derived model parameter 
of Teff. However, the fact that these uncertainties are ap- 
proximately reproduced provides evidence that the SVD 
method is valid. 

No identified mode is included with the binary system 
constraints (solid lines). The error box for the binary 
system does not shrink while reducing the errors in the 
radius (left panel). However, it shrinks when the error in 
the observed effective temperature is reduced, reproduc- 
ing accurately the observational ey^fj of (T(Teff) = 200, 
100, and 50 K. a(L*) does not decrease in either panel, 
because the mass is well-determined for the binary sys- 
tem and provides this strong constraint on L^,. 

In all cases the constraints provided by the binary sys- 
tem without an identified mode on the parameters of the 
pulsating star are more effective than those from the sin- 
gle star when an identified mode is included. 

5. INCORRECT MODE IDENTIFICATION 

Asteroseismology is only useful when the observed fre- 
quencies have been correctly identified with their mode 
geometry. Although techniques exist to do this, is there 
a way to test if this identification is correct? In this sec- 
tion we show that an incorrectly identified mode can be 
diagnosed if the pulsating star is in the binary system, 
but not as a single star. We subsequently show that we 
can in fact constrain and/or identify the mode correctly 
by studying the distribution of model solutions. 

We have shown that the oscillation mode is an im- 
portant observation to have (Sect. [3]) but including it 
as a constraint for the binary system did not improve 
the pulsating star's parameter uncertainties nor shrink 
the error box in the H-R diagram (Sect. |4]). This im- 
plied that the mode is somewhat redundant information 
for defining the global quantities, and we concluded that 
the observed oscillation frequency could possibly be used 
in a different way. The oscillation mode occupies the top 
rows of the U"^ matrix in the single star (Fig. |4|) and in 
the binary case, indicating that it is an observable that 
tightly constrains the stellar model. If the mode were 
incorrectly labeled, it should have a significant impact 
on the stellar parameter solution. We have also shown 
that the observational constraints from the binary sys- 
tem yield very well-determined parameters without using 
the oscillation mode (Figs. [71 [5]). Therefore, we can com- 
pare the model solutions when we use a set of observables 
with and without the oscillation mode as a constraint. 
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Fig. 8. — Maximum projections of the two-dimensional error ellipses defining the region of parameter space where the solution lies for 
the binary system, using OSBB-I-AS2 (black) and OSbs+ASS (grey), and the observational errors given in the third column of Table [21 




Fig. 9. — The theoretical error boxes for the model parameters of luminosity and effective temperature. The dashed/continuous lines 
represent the results for the single star (OSs+ASl)/binary system (OS^s). The left/right panels show the results while reducing the error 
in observed R/T^g. 



We perform simulations to test the effect of recovering 
the stellar parameters when we identify a mode incor- 
rectly. The model that we test has the true input param- 
eters (P) that are given in Table[T]with the corresponding 
model observables (B) from Tabled We generate a set of 
"real observations" {yi) by adding random errors to the 
true model observables: yi = Bi + eiri, where is a ran- 
dom number drawn from a Gaussian distribution with a 
mean of and standard deviation of 1 . To recover P from 
yi, we first make a guess of the solution Pg and then use 
equation ([2]) to recover the true parameters. We do this 
for four different scenarios for both the single star and 
the binary system, without including photometric infor- 
mation {OSg/EB and 0S5/£;s+ASl): (1) using the non- 
seismic measurements without a mode (black diamonds 



in Figure Uni OSs/eb), (2) using the non-seismic mea- 
surements and a correctly identified mode (red crosses 
in Figures [TOl and [TTl OSs/^b+ASI), (3) using the non- 
seismic measurements and an incorrectly identified mode 
— the wrong degree £ (blue crosses), and (4) using the 
non-seismic measurements and an incorrectly identified 
mode — the wrong radial order n (green crosses). 

The recovered parameters for 10,000 realizations are 
shown in Figures [TU] and [TT] We show only the four pa- 
rameters of the pulsating star discussed in Section 14.11 
The left/right panels show the results for the single star 
(OSs/OSs+ASl)/binary system (OSbb/OSeb+ASI). 
In each of the panels the black square shows the initial 
guess of the corresponding parameters, while the dotted 
lines show the model values from Table [T] 
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For the single star the left panels of Figure [TU] show 
that the non-seismic observables (OSs, black diamonds) 
do not constrain the parameters of the system. Once 
the mode is included (OS5+ASI), the solution becomes 
tightly constrained to a small range of parameter values. 
When the mode is correctly identified (red crosses) the 
recovered parameter range is correctly centered on the 
intersection of the lines. When the mode is incorrectly 
identified (blue or green crosses) the recovered parameter 
range is incorrect. It is possible to discard the green so- 
lution (the incorrect radial order) because the recovered 
values of X are outside of the typically accepted values. 
However, it is not possible to determine if the blue or 
the red values represent the true input parameters. This 
means that we are either heavily reliant on obtaining a 
correctly-labeled mode, or that the formal uncertainties 
are unrealistic as uncertainties because the systematic 
error is completely ignored. 

For the binary system the right panels of Figure [10] 
show that the parameter uncertainties agree with those 
shown in Figure [T] the l-cr uncertainty is shown by the 
error-bar. The three solutions with the mode all lead to 
well-constrained parameter ranges. The solution without 
the mode (black diamonds) can not be seen very clearly 
because the values coincide with the red solution, im- 
plying that it is the correct identification. The blue and 
the green solutions can be discarded simply because they 
are not in agreement with the black solution i.e. there is 
something incorrect about the mode-identification. 

To highlight the inconsistency in the solutions. Fig- 
ure [11] shows the fitted age for the three solutions in- 
cluding the mode (y-axis) versus the solution without 
the mode (x-axis). If the mode were correctly identified 
we would expect the solutions to approximately follow 
the bisector. This is precisely what we see in both pan- 
els. When identifying the mode with an incorrect de- 
gree (blue) or the incorrect radial order (green) there is 
a systematic offset, which implies an incorrectly labeled 
mode. The difference between the left and the right pan- 
els is that we included photometric information in the 
latter (OS£;s-|-AS2/AS3): here the offset is much clearer 
making it easier to detect the incorrectly-identified mode. 

This question can be alternatively posed to ask if we 
can use the distribution of model solutions to actually 
identify the oscillation mode. Figure [T^] shows a repre- 
sentative sample of the distribution of model solutions 
of similar simulations for the recovered age parameter. 
Each of the panels shows the recovered age without 
including the mode minus the recovered age using an 
arbitrarily-identified mode. We show just the results for 
the most ambiguous solutions, because all of the other 
mode identifications led to very clear discrepancies be- 
tween the two solutions. For the simulations our obser- 
vations come from a stellar model with a value of con- 
vective core overshoot parameter aov = 0.3. We show 
results when we assume that the inverting model has a 
value of Q!ov = 0.3 (correct) and 0.0 (incorrect). The 
top panel shows the results when the mode identification 
is correct {i,m,n) ~ (1,-1,9). The middle and lower 
panels show the results when the mode is labeled incor- 
rectly: (0,0,9) and (1,0,9). In both panels we can see 
that even with some errors in our assumptions about the 
model, we can safely constrain, or even label, the mode 
with the correct identification from the top panel. This 



docs not imply that the identied mode is insensitive to 
aov, simply that by comparing solutions from the same 
model any possible systematic error is eliminated. 

6. CONCLUSIONS 

Seismology of S Scuti stars has had limited success 
due to the inability to obtain a unique model of the 
star and more importantly the difficult task of mode- 
identification due to rotation, evolution and a deficiency 
in our understanding of mode-excitation mechanisms. A 
possible method for overcoming these obstacles, as sug- 
gested by many authors, is to study pulsating stars in 
detached eclipsing spectroscopic binary systems. In this 
paper we used singular value decomposition (SVD) to 
quantify the advantages of studying pulsating stars in 
eclipsing binaries over single field stars, by comparing 
the parameter uncertainties and the ability to detect an 
incorrectly-labeled mode in both cases. 

We have demonstrated that SVD is a powerful tool 
to assess the parameter constraints in various systems 
considering the available measurements. This method 
assumes that we are relatively close to the parameter 
solution so that the model derivatives can be described 
locally as linear. This implies that the results are not 
extendable to other values where we can no longer ex- 
trapolate linearly to estimate observables. We consid- 
ered the set of parameters shown in Table [1] and deviat- 
ing much from these values will require new calculations. 
This work also imposed various assumptions and they 
must be taken into account. For exa mple, the derivatives 
are c alculated from the ASTEC ((Christensen-DalsgaardI 
|2008[ ) stellar code, and they may vary (hopefully only 
slightly) from one code to the next. 

We have shown that as the measurement errors on the 
photometric information are varied, the roles of the ob- 
servables change especially for the single star. In par- 
ticular, R and (r — i) play similar roles for determining 
M, and (i - z) and [M/H] for determining Z. The dis- 
tance d is uniquely determined by tt when no photometric 
information is included, while the colors and magnitude 
contribute to determining d when these are included, pri- 
marily due to the error on the parallax. When no oscil- 
lation mode is included, then all of the observables are 
important for constraining the parameters even poorly, 
and the inclusion of a correctly identified mode replaces 
the role of R for determining M. 

For the binary system the photometric information 
does not have as much impact as it does for the sin- 
gle star, partially because these observations are blended 
values resulting from the two components. However on 
reaching precisions of 0.003 mag, then the photomet- 
ric observables begin to dominate the roles of Ra, TcsAi 
[M/H], and Tba for determining the stellar parameters. 
In particular, (r — i)EB yields similar information to 

TcffA- 

For a single star system correctly identifying an oscil- 
lation mode is necessary to reduce the parameter uncer- 
tainties to values that are useful. On the other hand, 
a binary system contains enough information, without 
having to use the identified mode to constrain the stellar 
parameters well. In fact, including the identified mode 
docs not improve the parameter uncertainties, emphasiz- 
ing its utility for learning something else about the star. 
This same trend can be seen for the error box in the H-R 
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Fig. 10. — Fitted parameters from simulations using OSeb and OSbs+ASI. The left/right panels show the results for the single 
star/binary. The color-code is as follows: black — no mode included; red, green and blue — including a mode, with various identifications: 
red — correctly identified; blue — incorrectly labeled degree; green — incorrectly labeled radial order. The intersection of the dotted lines 
are the true values of the input parameters, the black square shows the values of the initial guess, and the error bars on the right panel are 
the l-tr uncertainties for the binary system (see Fig.O. 
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Fig. 11. — Fitted age from the simulations when a mode is (y-axis) and is not (x-axis) included as an observational constraint. The color 
code is the same as Fig. 1101 and the left/right panels show the results for OSss+ASl/ASS. 
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Fig. 12. — Fitted age from simulations using OS+AS2 and OS+AS3. The left/right panels show the results for inversions using the 
correct /incorrect model (aov = 0.3/0.0). From top to bottom, the simulations assume a mode identification of (£,in,n) = (1,-1,9) 
(correct), (0,0,9) (incorrect) and (1,0,9) (incorrect). 
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diagram: it does not shrink for the binary system when 
the mode is included or the errors in Ra are improved. 
For the single star, however, improving R leads directly 
to shrinking the error box in the luminosity axis only 
when a correctly labeled mode is included. 

Other parameter values were also considered in this 
study; we varied the primary mass of the star (1.80 Mq 
and 2.50 Mq), the mass ratio of primary to secondary 
(1.10 and 1.60), the evolutionary stage of the primary 
star (central hydrogen mass fraction of 0.50 [MS] and 
0.15 [end of MS]), and the metallicity of both stars (0.020 
and 0.035). For a single star identifying a mode has more 
impact on the reference model than the more evolved 
model and for both the single star and the binary sys- 
tem the higher metallicity model benefited more from the 
photometric information. However, changing the size of 
the pri mary mass did n ot change the uncertainty depen- 
dencies (jCreevevll2008l ). 

Because mode-identification in 6 Scuti (and other) 
stars is a complicated task, we addressed whether an 
incorrectly-labeled mode could be correctly diagnosed 
using the observational constraints provided by the sin- 
gle star and the detached eclipsing binary system. We 
used simulations to test the effects of recovering stellar 
parameters with arbitrarily-identified oscillation modes. 
We found that for the single star, it is necessary to iden- 
tify the mode unambiguously, because the recovered so- 
lutions using an incorrect or a correct mode identifica- 
tion are not distinguishable. Meanwhile for the binary 
system, an incorrectly-identified mode can be diagnosed 
by simply comparing the parameter solutions when the 
mode is included and not. If the mode is correctly iden- 
tified, the solutions are in agreement. 

Taking this discussion slightly further, in the binary 
system we also showed that we could in fact tightly con- 
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strain (and even identify correctly) the pulsation mode 
parameters by comparing the distribution of model solu- 
tions with and without a mode. The model solutions are 
in agreement when the mode is correctly identified even 
if the assumptions on the interior physics arc slightly 
incorrect. This type of test and this method of mode- 
identification has not yet been done with available ob- 
servational data. 

It is clear that studying pulsations in binary systems 
has its complexities obscrvationally, however, the extra 
effort is definitely worth it. Apart from the advantage of 
using the screening effect of the pulsations during eclipse 
to help identify the mode, this study emphasizes the 
power that the binary measurements have for determin- 
ing fundamental parameters and subsequently constrain- 
ing and even identifiying the mode. 

In this paper we use the SVD technique to compare 
the constraints for single stars and binary systems. We 
note, however, that with the flood of data from space- 
based missions, this type of study can be subsequently 
extended to investigate which detached eclipsing binary 
systems provide the best astrophysical laboratories. 
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